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Abstract 


Compressive  measurement  in  speetroscopic  and  intcrfcrometrie  systems  has  been 
demonstrated.  Computational  spectroscopy  studies  resulted  in  the  development  of  the 
coded  aperture  snapshot  spectral  imager  (CASSI)  and  eoherenee  measurement  studies 
resulted  in  the  development  of  eompressive  holography.  We  developed  the  theoretical 
underpinnings  of  eoded  aperture  snapshot  speetral  imaging  and  designed  two  proof-of- 
coneept  systems:  a  pushbroom  system  (SmaeM)  and  a  snapshot  system  (MaeSim).  Iso¬ 
lated  growth  of  nanoseale  earbon  nanotube  (CNT)  islands  has  been  accomplished  as  well 
as  growth  via  a  Pt  eatalyst,  ineluding  microstructural  validation  of  the  CNT  structure  us¬ 
ing  transmission  electron  mieroseopy.  The  use  of  Pt  eatalyst  during  the  CNT  growth  ra¬ 
ther  than  the  more  traditional  CNT  eatalyst  materials  is  benefieial  due  to  the  extensive 
study  of  SERS  (Surface  Enhanced  Raman  Speetroseopy)  and  new  biologieal  application 
areas  are  enabled  due  to  Pt’s  approved  use  for  in-vivo  applications. 

2.  Objectives: 

The  initial  goals  of  this  program  were 

1.  To  design  and  demonstrate  a  platform  for  hypcrspcetral  Raman  mieroseopy 
based  on  results  from  previous  AFOSR  supported  work  in  computational  sensor 
design. 

2.  To  use  the  computational  hyperspectral  microscope  to  explore  spatial  and 
spectral  signatures  of  Raman  signals  generated  on  patterned  substrates  of 
earbon  nanotubes. 

3.  To  functionally  outline  optical  signals  and  nanotube-analyte  interactions 
capable  of  enabling  advanced  chemical  sensors  and 

4.  To  detail  charge  transport  and  molecular  interactions  in  CNT  deviecs  for 
sensing  and  photonic  devices. 

Supplemental  funding  was  added  to  the  program  in  2009  through  the  DARPA  MTO 
effort  “quantum  eoherenee  sensing.”  The  objective  of  the  extended  program  was  to 
explore  the  possibility  of  direet  measurement  of  optieal  eoherenee  functions  using 
quantum  interferometry.  Study  of  coherence  sensors  resulted  in  the  diseovery  of 
compressive  holography. 

3.  Status  of  effort: 

This  program  resulted  in  major  advances  in  the  use  of  compressive  sampling 
theory  to  reeonstruet  3D  data  volumes  from  snapshot  2D  measurements.  According  to 
previous  theory,  3D  reconstruction  required  temporal  scanning.  Under  this  program  we 
demonstrated  3D  speetral  data  cube  reconstruction  from  a  single  eoded  2D  measurement. 
Under  the  extended  effort  focusing  on  coherence  sensing,  we  further  demonstrated  3D 
tomographic  image  reconstruction  from  a  single  2D  Gabor  hologram.  These  develop¬ 
ments  hold  substantial  promise  for  applications  in  hyperspectral  remote  sensing  and  laser 
radar. 


A  Hadamard  Aperture  Microscopic  Multiple  Excitation  Raman  spectrometer  has 
been  designed,  eonstrueted,  and  its  operation  verified.  We  developed  the  theoretical  un¬ 
derpinnings  of  the  aperture  eoded  speetral  engine  and  designed  two  proof-of-eoneept  sys¬ 
tems:  a  pushbroom  system  (SmaeM)  and  a  snapshot  system  (MaeSim).  The  pushbroom 
system  is  based  on  a  statie  aperture  eoded  single  disperser  design  that  is  translated  in  or¬ 
der  to  capture  a  3D  eube  (2D  spatial,  ID  speetral).  The  snapshot  system,  however,  is  a 
dual  disperser  architecture  whieh  uses  advanced  signal  processing  techniques  (compres¬ 
sive  sensing)  to  measure  a  3D  eube  in  a  snapshot.  The  basic  snapshot  design  is  called  a 
Coded  Aperture  Snapshot  Speetral  Imager  (CASSI).  Imagers  based  on  this  design  com¬ 
putational  sensors  that  are  able  to  recover  a  3D  dataeube  representing  spatial  and  speetral 
information  at  every  pixel  in  an  image  of  a  scene  using  just  one  snapshot  2D  projection. 
CASSI  instruments  utilize  a  eoded  aperture  and  one  or  more  dispersive  elements  to  mod¬ 
ulate  the  optieal  field  from  a  seene.  A  deteetor  captures  a  two-dimensional,  multiplexed 
projection  of  the  three-dimensional  dataeube  representing  the  seene.  The  nature  of  the 
multiplexing  performed  depends  on  the  relative  position  of  the  eoded  aperture  and  the 
dispersive  elcment(s)  within  the  instruments.  Experimental  results  from  our  proof-of- 
eoneept  implementation  of  a  hyperspeetral  mieroseope  have  been  obtained  for  applica¬ 
tions  in  addition  to  simple  Raman  spectrometers. 

Isolated  growth  of  nanoseale  CNT  islands  has  been  accomplished  as  well  as 
growth  via  a  Pt  catalyst,  ineluding  mierostruetural  validation  of  the  CNT  structure  using 
transmission  electron  microscopy.  Plasma  enhanced  ehemieal  vapor  deposition  has  been 
utilized  to  grow  vertically  aligned  multiwall  carbon  nanotubes  (CNTs)  via  a  platinum 
catalyst.  A  thermally  oxidized  layer  of  SiC>2  served  as  a  diffusion  barrier  to  prevent  Pt 
silieide  formation  and  was  shown  to  aide  the  formation  of  well  dispersed  catalytic  nano- 
partieles.  Adjustment  of  the  eatalyst  film  thiekness  was  effective  for  varying  the  diame¬ 
ter  and  density  of  synthesized  CNTs.  TEM  analysis  revealed  a  bamboo-like  structure 
with  eatalyst  particles  in  the  tips,  root,  and  within  internal  compartments,  indicating  that 
the  eatalyst  was  in  a  liquid  state  during  growth.  The  use  of  Pt  eatalyst  during  the  CNT 
growth  rather  than  the  more  traditional  CNT  eatalyst  materials  is  beneficial  due  to  the 
extensive  study  of  SERS  (Surfaee  Enhanced  Raman  Spectroscopy)  utilizing  Pt.  In  addi¬ 
tion,  new  biological  application  areas  are  enabled  by  using  Pt  due  to  its  approved  use  for 
in-vivo  applications. 

This  program  was  extended  for  1  year  under  DARPA  MTO  seedling  support  focus¬ 
ing  on  coherence  measurement  and  compressive  sampling.  One  goal  of  this  portion  of  the 
study  was  to  analyze  the  feasibility  and  physical  implementation  of  quantum  coherent 
optieal  sensors.  As  part  of  the  seedling  effort,  we  performed  a  eareful  theoretical  study 
on  utilizing  quantum  coherence  to  perform  optieal  coherence  measurement,  and  identi¬ 
fied  a  physical  system  and  device  concept  to  realize  coherence  measurement  using  quan¬ 
tum  interactions.  The  debate  on  the  potential  advantage  of  utilizing  quantum  systems  for 
direct  measurement  of  optieal  coherence  between  two  fields  with  better  SNR  than  the 
classical  limit  has  not  quite  been  conclusively  settled.  Based  on  the  technologies  availa¬ 
ble  today,  physical  realization  of  such  coherence  measurement  remains  a  major  chal¬ 
lenge.  The  second  goal  of  this  portion  of  the  study  was  to  explore  the  use  of  elassieal  de¬ 
tectors  to  measure  compressive  projections  of  coherence  functions.  This  portion  of  the 
study  was  extremely  successful,  as  demonstrated  in  the  development  of  compressive  ho- 


lography.  We  demonstrated  independent  3D  tomographie  data  cube  of  various  samples 
using  compressive  holography. 


4.  Accomplishments/New  Findings: 


SmacM  Multispectral  Microscope 

A  second  generation  scanning  (pushbroom)  multispectral  aperture  coded  micro¬ 
scope  (SmacM)  applied  to  fluorescence  studies  has  been  constructed[l]. 


Base  tine 
Detector 


Flat 

Mirror 


Tube 

Relay  Lens 


20x-100x 

Ob)ecttve 


Intermediate 
Microscope 
tmage  Plane 


Intermediate 
Coupler 
tmage  Plane 


D 


Sample 


Scanning 

Stage 


Figure  1:  (Left)  Second-generation  system  schematic.  (Right)  Second-generation  system  photograph. 


An  aperture-  coded  spectrometer  is  interfaced  with  one  of  the  exit  ports  on  the 
Zeiss  Axioplan  2  microscope  and  translated  perpendicular  to  the  dispersion  direction  of 
the  spectrometer.  Coupling  optics  spatially  de-magnify  the  field  size  at  the  exit  port  of  the 
microscope  such  that  it  more  accurately  maps  to  the  maximum  field  size  accepted  by  the 
aperture  coded  spectrometer.  Translation  of  the  spectrometer  is  completed  using  a  three- 
axis  stage  controller  that  interfaces  with  a  Newport  ESP  300  motion  controller.  The  z- 
foeusing  capability  for  the  microscope  is  automated  using  another  three-axis  motion  con¬ 
troller.  The  Zeiss  Axioplan  2  microscope  was  modified  (sample  chamber  and  condensing 
optics  were  completely  removed)  for  laser  excitation  fluorescence  studies. 


Figure  2:  At  the  intermediate  image  plane  of  the  microscope  where  the  image  is  de-magmfied  wnth  a  CCTV 
0.35x  Nikon  coupler  and  interfaced  with  an  aperture -coded  spectrometer. 


Two  focal  planes  arc  used  in  the  setup  for  SmacM.  One  CCD  is  located  at  the 
baek-end  of  the  aperture  coded  spectrometer  and  the  baseline  camera  is  located  at  the 
second  output  port  of  the  microscope.  The  baseline  camera  is  a  SPOT  CCD  made  by  Di¬ 
agnostic  Instruments  is  a  multispeetral  camera  that  measures  three  spectral  channels  from 
an  image  using  a  liquid  erystal  tunable  filter.  Both  of  the  CCDs  can  be  seen  in  Figure  2. 


Bichromatic  Source 

As  a  proof-of-conccpt  test,  a  quartz  on  chrome  64  x  64  grid,  5.4  um  mask  with  the 
letters  “DISP”  is  illuminated  with  a  lOx  beam  expanded  Helium  Neon  and  green  Helium 
Neon  laser  profile.  Each  slice  in  the  data  cube  takes  less  than  a  second  to  collect.  A  total 
of  223  steps  are  collected  over  a  12mm  extent.  Figure  3  represents  slices  from  the  data 
cube  collected  from  a  spectrally  bichromatic  source. 


(a) 


D! 


(f) 


Figure  3:  Dataeube  generation  fora  two-lined  souree  illumination  of  a  chrome  on  quartz  mask  (D1SP).  (a) 
Image  at  spectral  channel  N  =  15.  (b)  Image  at  spectral  channel  N  =  97  showing  a  mask  pattern  obtained 
from  the  HeNe.  (c)  Image  at  speetral  channel  N  =  150.  (d)  Image  at  spectral  ehannel  N  =  248  showing  the 
full  mask  pattern  from  the  green  HeNe  (e)  Image  at  spectral  channel  N  =  350  (0  Baseline  2D  spatial  image 
obtained  with  the  SPOT  eamera  with  a  different  field  of  view. 


A  plot  of  a  row  of  voxels  at  a  particular  spatial  location  through  the  data  cube 
shows  two  delta  like  peaks  at  the  wavelengths  corresponding  to  the  helium  neon 
(632. 8nm)  and  green  helium  neon  peaks  (543. 5nm).  Speetral  slices  appear  in  the  green  at 
high  speetral  channels  and  red  at  lower  spectral  channels  for  the  reconstructed  data  cube. 


Figure  4;  Plot  for  a  row  of  voxels  in  a  data  cube.  Peaks  from  the  green  HeNe  laser  source  and  HeNe  laser 

source  are  shown. 


Broadband  Source 

Broadband  sources  were  also  tested  with  SmaeM.  Three  types  of  fluoreseent  mi- 
erospheres  (Fluosphcres)  from  Invitrogen  were  tested.  The  fluoreseent  mierospheres  were 
prepared  in  a  10:1  solution  with  de-ionized  water.  A  droplet  from  a  poly-l-lysine  solution 
is  plaeed  on  the  microscope  slide  before  adding  the  fluospheres  to  increase  the  adhesion 
process  of  the  fluospheres  to  the  microscope  slide.  De-ionized  water  is  then  used  to  wash 
the  droplet  off  of  the  microscope  slide  for  image  acquisition  with  the  instrument.  The  flu¬ 
ospheres  tested  consisted  of  0.04um  red-orange  fluoreseent  (565/580)  earboxylatc- 
modifed  mierospheres,  0.2um  orange  fluoreseent  (540/560)  earboxylate-modifed  miero¬ 
spheres,  and  1.0  m  nile  red  fluoreseent  (535/575)  earboxylate-modifed  mierospheres. 
Spectra  for  the  fluoreseent  mierospheres  were  first  collected  with  a  conventional  slit 
spectrometer,  Ocean  Optics,  as  a  baseline  test.  Results  with  the  aperture  coded  spectro¬ 
meter  were  then  compared  with  the  baseline  results. 


Ocean  Optics  spectrometer  Invitrogen  Fluosphere  spectra 
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Figure  5:  Overlay  of  fluorescent  microsphere  emission  spectra. 


Figure  6:  Data  cube  generation  for  0.2um  fluorescent  microsphcres  excited  by  532nm  laser  light,  (a)  Spec¬ 
tral  channel  N  =  268.  (b)  Spectral  channel  N  -  222  (c)  Spectral  channel  N  =  205.  (d)  Spectral  channel  N  = 
128.  (e)  Spectral  channel  N  =  10.  (f)  Baseline  image  (right)  (a)  Plot  along  the  spectral  channels  of  the  data 
cube  for  an  “on”  pixel,  (b)  Plot  along  the  spectral  axis  of  an  “off”  pixel. 


Figure  6  shows  the  results  obtained  from  the  0.2um  fluorescent  mierosphere  test.  A  base¬ 
line  image  is  used  to  compare  the  spatial  resolution  of  the  reconstructed  spectral  slices  for 
a  given  data  set.  When  plotting  a  row  of  voxels  at  a  spatial  location  defined  by  the  fluo- 
sphere,  a  fluorescence  spectral  profile  is  obtained.  This  spectral  profile  can  be  compared 
to  that  obtained  from  a  slit  spectrometer  as  seen  in  Figure  5. 


Spatial  and  Spectral  Reconstructions 

Testing  spatial  and  spectral  reconstructions  from  fluorescent  mierospheres  is  a 
baseline  test  conducted  with  eonfoeal  imaging  systems.  The  next  step  in  verifying  the  ap¬ 
plicability  of  the  SmacM  system  to  fluorescence  microscopy  studies  is  to  analyze  fixcd- 
eellular  assays.  Data  with  SmaeM  was  taken  from  a  FluoCell  microscope  slide  #2  from 
Invitrogen.  The  slide  contained  bovine  pulmonary  artery  endothelial  cells  stained  with 
red-fluoreseent  Texas  Red-X  phalloidin  for  labeling  F-aetin,  mouse  monoclonal  anti¬ 
alpha-tubulin,  green-fluorescent  BODIPY  FL  goat  anti-mouse  IgG  was  used  to  label  mi¬ 
crotubules  and  blue-fluoreseent  DAPI  was  used  to  label  the  nuclei.  Due  to  the  sensitivity 
of  SmaeM  only  two  (red  &  green)  out  of  the  three  fluorescent  signals  could  be  detected. 
The  slide  was  again  excited  with  532nm  light  and  imaged  through  a  50x  objective,  NA  = 
0.5  and  a  0.35x  coupling  optic. 


Figure  7:  Invitrogen  cell  slide:  Emission  spectra  from  BODIPY  and  IgG  fluorescent  labels  in  the  bovine 

pulmonary  artery  endothelial  cell. 


Figure  7  describes  the  baseline  emission  profiles  from  two  of  the  fluorescence  signals 
from  the  Invitrogen  slide. 
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Figure  8:  Data  cube  generation  for  the  bovine  cell  on  a  microscope  slide  excited  by  a  532nm  laser  source 
(a)  Spectral  channel  N  =  310  (b)  Spectral  channel  N  =  275.  (c)  Spectral  channel  N  =  228.  (d)  Spectral 
channel  N  =  121.  (e)  Spectral  channel  N  =  10.  (0  Baseline  image  for  the  bovine  endothelial  cell  with  the 

SPOT  camera. 


Figure  8  describes  the  reconstructed  data  cube  obtained  for  the  bovine  endothelial 
cell.  Compared  to  the  baseline  image  the  spatial  resolution  of  the  reconstructed  slices 
eould  be  improved.  The  two  fluorescent  signals  are  distinguished  within  the  data  cube. 


Target  Classification 

Under  a  previous  AFOSR  grant  #3136057,  wc  have  developed  a  spectrally  adap¬ 
tive  signal  segmentation  algorithm  that  allows  the  direct  classification  of  targcts[2].  To 
demonstrate  the  classification  capabilities  of  the  algorithm  we  have  used  as  a  test  bed  a 
DD-CASSI  imager  incorporated  into  a  Zeiss  AxioObscrvcr  A1  fluorescence,  microscope 
(MacSim),  see  Figure  9 
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Figure  9:  (A)  System  layout  with  a  Zeiss  AxioObserver  Al  upright  microscope.  (A  1 )  Zemax  layout  for  the 
Arm  1  in  the  DD-CASSI  optical  design. (B)  Experimental  setup.  (Bl)  DD-CASSI  interfaced  with  the  output 

port  of  an  upright  microscope. 


In  collecting  the  data  a  spectrum  is  reconstructed  at  every  pixel  location  in  the  im¬ 
age  using  a  total  variation  (TV)  minimization  compressive  sensing  algorithm.  The  object 
being  imaged  contains  polystyrene  microsphcrcs  with  ten  unique  spectral  emissions 
shown  in 

Figure  10.  A  total  variation  (Tlr)  minimization  algorithm  is  then  adapted  for  im¬ 
age  segmentation  of  the  recorded  spectral  image.  To  segment  the  image,  a  priori  know¬ 
ledge  of  the  spectral  signatures  expected  in  the  image  or  used  during  the  segmentation 
process.  The  instrument  in  Figure  9  records  3 1  spectral  channels  spanning  the  spectral 
range  between  450nm  and  730nm  with  lOnm  resolution. 
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Figure  10:  Spectral  profiles  for  different  fluorescent  microspheres 


A  two-step  iterative,  shrinkage  and  thresholding  algorithm  (TwlST)  has  been 
adapted  to  estimate  the  3D  (x,y,A.)  data  cube  from  a  2D  snapshot  taken  by  the  CASSI  sys¬ 
tem.  To  identify  the  desired  signatures  a  spectrally  adaptive  image  segmentation  algo¬ 
rithm  based  a  total  variation  minimization  compressive  sensing  algorithm  was  applied. 
The  result  is  a  3D  data  eube,/fjc,v,  X),  estimate  projected  into  a  lower  dimensional  tex¬ 
ture-based  data  cube,  0(x,y,n).  Each  texture  is  associated  with  a  spectral  signature,  n, 
unique  to  the  sample.  The  new  spectrally  textured  data  cube  records  the  type  (n)  of  each 
pixel  registered,  as  well  as,  the  brightness  value  associated  with  pixel  location.  In  Figure 
1 1 ,  a  transformation  from  f  to  0  is  described.  Each  n  value  corresponds  to  an  associated 
spectral  signature  for  the  cube  layer. 
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Figure  1 1 :  Transformation  of  the  f(x,y,  a)  data  cube  into  a  lower  dimensional  0(x,y,n)  datacubc 

The  transformation  of  f  to  0  is  represented  by 

f=  W  0 

O'  =  argmin||g-  HW6|^  +  rd>n,($) 

a  2 

<M«> = =I|«  v">-.  J 
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Again  a  total  variation  minimization  approach  is  used  to  converge  to  an  image 
with  the  sparsest  gradient.  In  the  above  equations,  H  is  a  positive,  semi-definite  matrix,  g 
is  the  intensity  recorded  at  the  focal  plane,  0  is  an  estimate,  t  is  the  regularizing  parameter 
that  minimizes  the  ill-posed  nature  of  the  problem,  and  tp  is  the  regularizing  function  that 
applies  a  gradient  along  the  rows  and  columns  of  the  cube  estimate.  Figure  12  provides 
an  intuitive  example  of  how  a  gradient  impacts  an  image.  Four  squares  are  used  to  model 
fluorescent  mierospheres  of  different  spectral  emission  and  different  intensities.  A  gra¬ 
dient  is  similar  to  an  edge  detection  filter  as  seen  in  the  diagram  below. 


Figure  12:  Demonstration  of  the  use  of  total  variation  minimization  to  produce  the  sparest  gradient  of  an 
image.  Flerc  four  squares  of  unique  color  are  reduced  to  outlines  of  the  squares. 

An  algorithm  flow  diagram  is  shown  in  Figure  13.  An  image  is  recorded  by 
CASSI  instrument  in  Figure  9  and  intensity  normalized.  The  1200x  1600  image  is 
cropped  to  a  local  region  corresponding  to  the  field  of  view  of  single  bead,  g.  For  the 
adaptation  of  TwlST  for  image  segmentation,  prior  information  in  the  form  of  a  norma- 


lized  spectral  database,  W,  and  a  cube,  C,  corresponding  to  a  spectral  impulse  response 
every  5nm  between  450nm-730nm  is  provided  to  the  algorithm. 
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Figure  13:  Image  segmentation  algorithm  flow  diagram 


Filtered  light  from  a  Xenon  arc  lamp  attached  to  the  microscope  with  a  bandwidth  from 
377nm  to  396nm  is  used  to  excite  the  polystyrene  microsphcrcs  that  make  up  the  sam¬ 
ples.  An  intermediate  image  plane  is  formed  at  the  output  ports  of  the  microscope.  The 
DD-CASSI  is  interfaced  to  the  microscope  via  the  output  port  and  captures  the  emission 
bandwidth  from  420nm-800nm.  A  baseline,  RGB  detector  is  placed  at  the  second  output 
port  of  the  microscope.  A  high  frame  rate  camera  that  allows  dynamic  scene  analysis  and 
a  cooled,  long  exposure  capable  camera  for  static  samples  can  be  used  to  record  the  out¬ 
put  of  the  CASSI. 

A  fluorescent  scene  recorded  by  the  CASSI  system  is  shown  in  Figure  914.  The 
image  contains  all  ten  bead  varieties. 


Figure  14:  Raw  DD-CASSI  image 


An  RGB  baseline  image  is  obtained  from  the  seeond  port 


Figure  15:  RGB  image  of  a  fluorescence  scene. 


It  is  a  challenge  to  distinguish  between  the  ten  different  bead  varieties  shown  in 
Figure  154.  Using  the  image  segmentation  algorithm  we  have  developed,  a  bead  identity 
map  (or  spectral  texture-based  map)  clearly  allows  one  to  distinguish  all  ten  bead  varie¬ 
ties,  as  shown  in  Figure  16. 


Figure  16.  Image  segmentation  results 

Figure  1 7,  displays  a  pseudo-color  estimate  from  reconstructed  data. 


Figure  17:  Pseudo-color  image  of  fluoresecnl  beads. 


In  summary,  coded  aperture  snapshot  spectral  imager  has  been  constructed  and 
applied  to  fluorescence  microscopy  applications.  A  total  variation  minimization  algo¬ 
rithm  has  been  applied  to  the  CASSI  snapshot  spectral  imager  to  allow  a  spectrum  to  be 
reconstructed  at  every  pixel  location  in  an  image.  The  algorithm  developed  allows  the  use 
of  image  segmentation  based  on  texture  to  be  used  for  fluorescence  bead  analysis.  Both 
static  and  dynamic  fluorescent  beads  have  been  analyzed. 
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Figure  18:  Plol  along  lhe  speelral  channels  for  a  given  pixel  loealion 


Carbon  Nanotube  Growth  on  Pt 

Sample  preparation  involved  cleaving  silicon  (100)  wafers  to  desired  size  and  ul- 
trasonically  cleaning  in  solvent.  The  substrates  were  then  dried  in  nitrogen  and  oxidized 
in  a  tube  furnace  at  1000°C  in  dry  02  for  2  hours  and  allowed  to  cool  overnight.  The  re¬ 
sulting  oxide  thickness  was  ~  200  nm.  Pt  thin  films  of  three  different  thicknesses  (15,  34, 
and  50  A)  were  deposited  on  the  oxidized  substrates  using  a  CHA  Electron-beam  Evapo¬ 
ration  system.  Pt-eoated  substrates  were  then  placed  in  a  915  MHz  M-PECVD  reactor 
and  pumped  down  to  <  10  mTorr  in  preparation  for  CNT  growth.  The  use  of  ammonia  or 
hydrogen  plasma  for  pretreatment  resulted  in  similar  CNT  growth  morphology  with  the 
exception  that  ammonia  pretreatment  yielded  slightly  better  vertical  alignment.  The 
presence  of  nitrogen  species  is  known  to  enhance  nueleation  during  MPECVD  growth 
creating  a  higher  density  of  CNTs  which  have  a  greater  interaction  via  van  der  Waals 
forces  promoting  growth  along  the  same  direction.  This  effect  would  be  enhanced  as 
CNTs  grow  longer.  Nevertheless,  for  the  remainder  of  this  report  we  will  simply  refer  to 
hydrogen  or  ammonia  similarly  as  reducing  gas  or  plasma. 


Prior  to  growth,  the  substrates  were  heated  to  1000°C  in  100-150  seem  of  reduc¬ 
ing  gas,  followed  by  striking  and  stabilizing  a  reducing  plasma  at  21  Torr  and  21  kW  of 
magnetron  input  power.  It  should  be  noted  that  the  growth  temperature  specified  is  the 
temperature  setting  of  the  substrate  heater  and  does  not  accurately  reflect  the  temperature 
of  the  sample  due  to  its  distance  from  the  thermocouple  measurement  device  and  also 
heating  induced  by  the  plasma  itself.  Substrates  were  then  pretreated  for  1  minute  under 
the  stabilized  reducing  plasma.  The  total  heat-up  and  pretreatment  time  was  not  more 
than  5  minutes  during  which  the  Pt  film  dewetted  into  discrete  nano-sized  islands.  Fol¬ 
lowing  pretreatment,  the  reducing  gas  flow  was  halted  and  a  mixture  of  methane  and 
ammonia  gas  with  a  gas  ratio  of  3  -  5  CH4/NH3  was  introduced  into  the  plasma  to  initiate 
the  growth  of  CNTs.  Total  flow  rates  did  not  exceed  200  seem  and  growth  times  did  not 
exceed  6  minutes. 

The  monitoring  of  nanopartiele  evolution  was  critical  for  the  successful  growth  of 
Pt-CNTs  as  the  conditional  process  space  was  found  to  be  much  narrower  than  that  for 
Fe-CNT  growth  with  similar  substrate  material  in  our  reactor.  For  instance,  Pt-CNT 
growth  did  not  occur  at  <  1000°C  no  matter  the  change  in  various  other  experimental  pa¬ 
rameters  including  growth  gas  ratio,  pretreatment  or  growth  time,  catalyst  thickness,  etc. 
Figure  19  presents  cross-section  SEM  micrographs  of  nanopartiele  evolution  following 
different  pre-growth  conditions. 


Figure  19:  Cross-section  SEM  images  showing  nanopartiele  evolution  after  reducing  plasma  pretreatments 
at  lOOO’C  for,  (a)  50  A  Pt  /  Si  after  1  minute,  and  fb)  50  A  Pt  /  2000  A  Si02  /  Si  after  1,  (e)  3  ,  and  (d)  10 

minutes. 


It  was  observed  that  without  the  SiO;  layer,  a  3-dimensional  network  of  discrete 
and  spherical  Pt  nano-islands  did  not  form.  When  the  oxidized  layer  was  present,  dis¬ 
crete  hemispherical  shaped  Pt  nano-islands  could  be  obtained  under  the  reducing  plasma 
although  an  increase  in  pretreatment  time  resulted  in  flatter  and  more  round  nano-islands 
Interestingly,  the  Pt  nano-islands  were  also  partially  embedded  in  the  SiO:  layer  and  the 
depth  of  partial  submergence  increased  with  pretreatment  time.  This  change  in  catalyst 
island  morphology  with  increasing  pretreatment  time  is  believed  to  have  resulted  in  in¬ 
creased  carbon  encapsulation  of  Pt  nano-islands  yielding  very  short  CNT  lengths  (<  200 
nm)  and  eventually  no  CNT  growlh  at  all  during  the  growth  stage.  However,  it  is  be¬ 
lieved  that  the  porous  silica  layer  may  help  isolate  and  prevent  agglomeration  of  the  Pt 
nano-islands  as  well  enhance  catalytic  efficiency  of  Pt  particles  and  thereby  promote  the 
growth  of  CNTs. 


X-ray  diffraction  (XRD)  has  previously  been  effectively  utilized  to  identify  ultra- 
thin  platinum  silieide  films  deposited  on  (100)  silicon  and  was  thus  used  to  investigate 
the  possibility  of  substrate-catalyst  interactions  in  this  study.  Phase  analysis  experiments 
were  conducted  via  a  Philips  X’Pcrt  PRO  MRD  HR  XRD  System  to  monitor  the  atomie 
bonding  configurations  at  the  interface  of  the  deposited  Pt  film  on  the  substrate  materials 
before  and  after  pretreatments.  The  X-ray  source  delivered  primarily  Cu  Ka  X-rays  with 
an  incident  wavelength  of  1.5405  A.  Weaker  diffraction  peaks  from  Cu  K.p  radiation 
were  also  observed  in  the  resulting  spectra  when  the  Cu  Ka  diffraction  signal  was  in¬ 
tense,  which  occurred  when  probing  the  Si  wafer.  As-deposited  Pt/Si  substrate  exhibited 
a  broad  peak  at  39.8°  which  corresponds  to  a  disordered  Pt  (1 1 1)  lattice  structure.  Inte¬ 
restingly,  after  a  pretreatment  time  of  one  minute,  the  Pt  (1 1 1)  peak  is  significantly  dimi¬ 
nished  in  the  XRD  spectrum  and  becomes  absent  after  10  minutes.  Simultaneously, 
many  small  peaks  attributed  to  various  crystallographic  orientations  of  platinum  silieide 
(JCPDF  #7 10523)  arise. 

Alternatively,  XRD  confirmed  that  the  SiO;  interlayer  acted  as  a  diffusion  barrier 
to  silieide  formation  during  the  pretreatment  conditions.  A  sharp  Pt  ( 1 1 1 )  peak  was  de¬ 
tected  after  a  1  min  pretreatment  in  reducing  plasma  indicating  that  the  dewetted  Pt  metal 
was  able  to  crystallize  and  form  the  hemispherical  nanopartieles  shown  in  Figure  19. 
However,  as  the  pretreatment  time  was  extended,  the  amount  of  crystalline  Pt  was  conti¬ 
nually  redueed,  and  after  10  minutes,  the  amount  of  Pt  (1 1 1)  phase  was  below  the  instru¬ 
ment’s  detection  limit.  No  other  peaks  besides  those  attributed  to  Pt  ( 1 1 1 )  and  the  Si 
substrate  were  detectable  within  the  acquired  range  contrary  to  what  was  observed  after 
pretreatment  when  the  SiCL  interlayer  was  absent.  The  loss  of  crystalline  Pt  at  longer 
pretreatment  times  may  be  due  to  sublimation,  or  diffusion  of  Pt  atoms  within  the  SiO; 
interlayer  creating  a  diffuse  impurity  within  the  SiC>2  lattiee.  The  submergence  behavior 
previously  discussed  and  depicted  in  Figure  19  favors  the  latter  hypothesis.  Ultimately, 
the  XRD  results  of  pretreated  CNT  substrates  reflect  observations  noted  in  the  SHM  im¬ 
ages  and  further  study  of  substrate-eatalyst  interactions  may  have  implications  for  pro¬ 
moting  the  efficiency  of  Pt-CNT  growth. 

Scanning  and  Transmission  Electron  Microscopy  of  CNTs 

High-resolution  cross-section  and  plan-view  SEM  images  of  successful  Pt-CNT 
growth  using  short  duration  (<  1  min.)  pretreatments  for  three  different  eatalyst  thick¬ 
nesses  are  shown  in  Figure  20. 

Catalyst  thickness  was  an  effective  parameter  for  controlling  density  and  CNT  di¬ 
ameter.  Given  the  relatively  high  surface  tension  of  Pt  compared  to  silicon  dioxide,  the 
formation  of  nano-islands  on  SiCL  is  expected  to  adhere  to  the  Volmer-Wcbcr  growth 
mode.  Therefore  as  a  result  of  pretreatments,  thinner  catalyst  films  resulted  in  a  higher 
density  of  smaller  nano-islands  and  thus  forming  a  dense  forest  of  Pt-CNTs  with  smaller 
diameters  (a,  b).  Likewise,  thicker  eatalyst  films  resulted  in  a  lower  density  of  larger  na¬ 
no-islands  which  formed  by  ripening  and  eoaleseenee.  These  formed  larger  diameter  Pt- 
CNTs  and  a  lower  density  film  (e,  f).  A  Pt  film  thickness  of  ~  8  nm  did  not  produce  suc¬ 
cessful  growth  of  CNTs  as  resulting  nanopartieles  were  too  large  to  support  CNT  growth. 
The  average  diameter  of  CNTs  synthesized  from  15,  34,  and  50  A  thick  Pt  catalyst  films 
were  14,  28,  60  nm  respectively.  In  all  eases,  Pt-CNTs  had  uniform  diameter  distribu¬ 
tions  which  did  not  exceed  a  range  of  ~  30  nm. 


Figure  20:  Cross-section  (left  column)  and  plan-view  (right  column)  SEM  images  of  Pt-CNT  growth  given 
the  follow  ing  Pt  film  thicknesses,  (a)  and  (b),  15  A;  (c)  and  (d),  34  A;  and  (e)  and  (0,  50  A 

SEM  images  also  clearly  show  catalyst  particles  in  the  tips  suggesting  a  tip  growth 
mode,  although  particles  were  also  seen  at  the  root  and  within  the  Pt-CNTs  (which  is  dis¬ 
cussed  further  in  the  TEM  analysis  section).  The  tip-growth  mode  has  recently  been 
shown  to  correlate  with  greatly  reduced  growth  rates.  This  is  associated  with  significant 
reshaping  of  catalyst  nanoparticles  and  premature  earbon  encapsulation  during  CNT  for¬ 
mation  due  to  a  weak  support-catalyst  interaction  and  fast  sclf-diffusivity  of  the  catalyst 
metal.  Encapsulated  Pt  particles  in  our  CNTs  were  typically  elongated  and  greatly  re¬ 
shaped.  Recently,  a  transition  from  tip-growlh  mode  to  a  base-growth  mode  for  random¬ 
ly  oriented  Pt-CNTs  has  been  reported  in  which  the  carbon  yield  was  significantly  im¬ 
proved.  Experimental  parameters  of  our  PECVD  system  were  controllably  varied  in  a 
design  of  experiments  in  an  attempt  to  extend  the  length  or  thickness  of  the  Pt-CNT 
films,  however  this  effort  was  unsuccessful.  Goals  of  future  work  include  modification 
of  substrate  materials  to  increase  catalytic  efficiency  and  longevity  to  sustain  growth  of 
longer  aligned  Pt-CNTs. 

Low  resolution  TEM  was  utilized  to  inspect  the  internal  structure  of  the  Pt-CNTs. 
Preparation  of  TEM  samples  consisted  of  sonicating  as-deposited  samples  in  a  1:1  mix¬ 
ture  of  20  mL  of  acetone  and  ethanol  for  20  minutes  to  remove  Pt-CNTs  from  the  sub- 


stratc  and  reduce  agglomeration.  Then  using  a  pipette,  drops  of  the  resulting  solution 
were  placed  on  a  carbon  holey  film  supported  by  a  copper  mesh  TEM  grid.  The  image  in 
Figure  21  clearly  shows  a  bamboo-like  structure  given  the  scries  of  hollow  and  periodi¬ 
cally  stacked  cone-shaped  compartments.  Again,  catalyst  particles  arc  easily  seen  at  the 
tips  of  the  CNTs  suggesting  a  tip-grow  th  mode.  However,  upon  further  inspection,  cata¬ 
lyst  particles  are  also  commonly  seen  at  the  root  of  the  CNTs  and  the  thickness  of  the  in¬ 
ner  walls  decrease  from  tip  to  root  along  the  hollow  compartments. 

The  orientation  of  the  cone-shaped  catalyst  particle  encapsulated  at  the  tip  of  the 
CNT  is  opposite  the  orientation  of  most  of  the  cone-shaped  compartments  along  the 
length  of  the  CNT  which  have  the  same  orientation  as  the  catalyst  particle  at  the  root.  In 
addition,  some  catalyst  is  occasionally  observed  within  compartments  along  the  length  of 
the  tube.  Furthermore,  the  catalyst  particles  encapsulated  at  the  tips  arc  significantly 
larger  than  the  diameter  of  the  corresponding  nanotubc.  This  suggests  that  during  high 
temperature  growth  the  catalyst  metal  is  in  a  liquid  state  and  increasingly  more  parts  of 
the  particle  will  be  sucked  into  and  toward  the  growth  front  of  nanotube  due  to  capillary 
forces.  Overall,  these  observations  for  the  present  Pt-CNT  growth,  suggest  that  the  cata¬ 
lyst  particles  anchored  at  the  root  by  the  oxidized  silicon  substrate  determines  the  growth 
of  the  nanotubc  via  a  base-growth  mode,  and  the  catalyst  particle  at  the  tip  is  actually  not 
active.  In  fact,  the  depletion  of  the  initial  platinum  particle  at  the  root  being  sucked  to¬ 
ward  the  tip  may  contribute  to  the  cessation  of  Pt-CNT  growth.  However,  more  characte¬ 
rization  is  needed  to  confirm  the  details  of  the  growth  mechanism. 


Figure  21 ;  ThM  images  exhibiting  bamboo-like  structure  consisting  of  hollow,  truncated  conical  com¬ 
partments  and  tip  of  Pt-CNT  with  encapsulated  Pt  particle 

Nanoscale  Carbon  Nanotube  Island  Growth 

Using  a  FIB  (Focused  Ion  Beam)  to  define  the  Pt  catalyst  area,  nanoscale  CNT  Isl¬ 
ands  can  be  grown  (Aknowledgement:  Anuj  Dhawan).  The  sample  shown  in  Figure  22 
was  grown  with  50  angstrom  of  platinum  and  a  150:50  Methanc:Ammonia  ratio  at  900C. 
These  islands  will  enable  controlled  sensing  measurements  with  well  defined  CNT  vo¬ 
lumes.  Future  work  will  utilize  Raman  spectroscopy  to  examine  these  Pt  catalyzed  CNTs 


with  varying  nano-island  dimensions  and  analytes. 


Quantum  Coherence  Sensor 

The  goal  of  this  portion  of  the  study  was  to  analyze  the  feasibility  and  physical  im¬ 
plementation  of  quantum  coherent  optical  sensors.  As  part  of  the  seedling  effort,  we  per¬ 
formed  a  careful  theoretical  study  on  utilizing  quantum  coherence  to  perform  optical  co¬ 
herence  measurement,  and  identified  a  physical  system  and  device  concept  to  realize  co¬ 
herence  measurement  using  quantum  interactions.  The  debate  on  the  potential  advantage 
of  utilizing  quantum  systems  for  direct  measurement  of  optical  coherence  between  two 
fields  with  better  SNR  than  the  classical  limit  has  not  quite  been  conclusively  settled. 
Based  on  the  technologies  available  today,  physical  realization  of  such  coherence  mea¬ 
surement  remains  a  major  challenge. 

Classical  interferometry  for  measuring  the  mutual  coherence  function  Ij2  is  illu¬ 
strated  in  Fig.  23a.  The  two  optical  fields  are  described  by  a  set  of  annihilation  operators 
bk0  and  bk]  w  ith  a  quasi-monoehromatie  spectrum  featuring  a  center  frequency  of  a )c, 
bandwidth  of  Ao),  and  relative  phase  shift  of  cp.  Wc  assume  that  the  spectrum  of  the  two 
fields  is  identical,  so  that  Ij,  =  T22  =  /„  and  (blb^-j  =  where  nk  is  the  mean  photon 
number  in  the  mode  k.  The  normalized  mutual  coherence  ^=Y j2//0  characterizes  the  level 
of  coherence  between  the  modes  bk0  and  bkr  The  irradianee  on  the  detectors  arc 

/l2M  =  /n[l±Re(^)]. 

Measurement  of  Ii2((p)  at  three  distinct  values  of  (p  enables  independent  estimation 
of  the  phase  and  amplitude  of  Ij2.  When  |£f  is  small,  the  variance  of  each  measurement  of 
/,  2  is  approximately  equal  to  the  variance  of  measurements  of  /0,  given  by 

cr,  =2rjI()(l  +  T]n)h(ojTA, 

where  t]  is  the  quantum  efficiency  of  the  detector,  n  =  l/(eh"'  —  1)  is  the  mean  occupa¬ 

tion  number  at  (oc  for  a  source  with  effective  temperature  0,  T  is  the  integration  time,  and 
A  =  nWo2  is  the  cross-sectional  area  of  the  optical  fields  assumed  to  be  in  Gaussian  mod¬ 
es  with  beam  waist  Wo.  The  “noise  equivalent  coherence”  corresponding  to  |ljZmin|2  =  of 

is  a  measure  of  the  minimum  detectable  mutual  coherence  in  the  classical  system.  |rj,niin|’ 


is  independent  of  and  the  signal-to-noise  ratio  (SNR)  in  estimating  rj,  degrades  as  £ 
decreases. 


Conventional  Photodetectors 


Quantum  Coherence 
(b)  Detector  in  Optical  Cavity 


Figure  23:  Schematic  for  sensing  mutual  coherence  using  (a)  a  beam  splitter  and  classical  irradiancc  dctcc- 

tors  and  (b)  quantum  coherence  detector. 

Our  model  quantum  eoherenee  detector  replaces  the  beamsplitter  and  the  two  pho- 
todeteetors  (Fig.  23b)  with  a  quantum  system  directly  coupled  to  the  two  optical  fields. 

As  a  physical  implementation,  we  consider  a  collection  of  three-level  atoms  in  a  small 
optical  cavity  whose  level  structure  is  shown  schematically  in  Fig.  24a.  This  system 
choice  provides  a  pair  of  stable  ground  states  (e.g.,  hyperftne  levels  of  an  atom  with  non¬ 
zero  nuclear  spin)  that  feature  long  eoherenee  times  and  negligible  inhomogeneous  broa¬ 
dening.  The  excited  state  (denoted  by  |2)  )  is  coupled  to  two  ground  states  (|0)  and  |1 )  ) 
via  Raman  transition.  The  transitions  between  the  two  ground  states  and  the  excited  state 
are  ehosen  to  require  orthogonal  polarization  corresponding  to  the  optieal  fields.  We  also 
consider  two  coherent  eontrol  beams  of  frequency  cut  and  a>z  described  by  the  operators 
dx  and  d2  interacting  with  the  same  atoms,  where  a  complete  eontrol  of  the  mutual  eohe¬ 
renee  is  available  by  manipulating  the  intensities  and  relative  phase  of  the  beams.  The 
cavity  is  necessary  to  enhance  the  coupling  between  the  optieal  fields  and  the  atom,  so 
that  the  coherent  interaction  necessary  for  the  eoherenee  of  the  field  to  drive  the  atomic 
transition  occurs  within  the  deeoherenee  time  of  the  atoms  due  to  spontaneous  emission. 


Figure  24:  (a)  Level  structure  of  the  atom  under  consideration.  The  lwo  ground  stales  |0)  and  |1 )  are 
coupled  via  an  excited  state  |2)  through  optical  dipole  transitions,  (b)  Bloch  sphere  representation  of  the 
atomic  syslem.  The  interaction  with  the  field  rotates  the  initial  state  around  the  axis  n,  located  on  the  x  -  y 
plane  making  an  angle  <p  with  the  x-axis,  by  the  angle  2|i7o|  T. 


The  atomie  population  dynamics  ean  be  modeled  using  standard  Jaynes-Cummings 
interaction  Hamiltonian  under  adequate  approximations.  Under  the  assumption  of  large 
detuning  ( Au  Ak  »  Aco)  and  considering  time  evolution  over  time  seales  much  longer  than 
the  optieal  bandwidth  (t  »  Aco),  the  probability  amplitudes  for  the  two  ground  states  (C|(/) 
and  C2(/))satisfy  a  simple  pair  of  equations 

C’o(t)  —  +  t!2(f)C  j(<). 

c\(t)  -  tfnooM  +  id.c.io, 


where  Si  denotes  the  AC  Stark  shifts  of  eaeh  level  due  to  interaction  with  the  optieal 
fields  and 


ft(0  =  £  9*029^  12 


is  the  (eomplex)  Rabi  frequency  of  the  system.  Under  large  detuning  assumption,  the  Ra- 
bi  frequency  ensemble-averaged  over  the  bandwidth  of  the  optical  fields  beeomes  time- 
independent  and  reduees  to 

Z0d02dl2  u>02“J\2  , 


Oo  =  ((!(<))  a 


■r,j(T). 


hr  A  u>ouJ\  t r 

where  da  denote  the  dipole  moment  between  the  two  atomic  states  |z)  and  |2)  ,  A  =  {  At) 
and  ?oi  =  ( coki )  arc  the  ensemble  average  values  of  the  detuning  and  the  optical  frequency, 


respectively,  taken  over  the  optical  field,  and  the  phase  difference  is  taken  to  be  cpk  =<vtlr 
.  Zo  is  the  impedance  of  vacuum,  and  ^is  the  finesse  of  the  optical  cavity.  We  further  ex¬ 
pressed  the  electric  field  using  the  field  annihilation  operators  and  used  the  definition  for 
rj2(r).  We  see  that  the  ensemble-averaged  Rabi  frequency  is  proportional  to  the  mutual 
coherence  fj2(r)  of  the  two  optical  fields.  If  the  atoms  arc  prepared  in  the  initial  state 


|¥(t  =  0)}  =  coso  1 0)  +  c^sinall). 


one  can  find  an  exact  solution  to  the  equations  that  describe  the  atomie  population  dy¬ 
namics.  When  the  two  levels  |0)  and  1 1 )  are  exactly  degenerate  in  the  presence  of  AC 
Stark  shifts  (So  =  Si),  the  system  undergoes  a  coherent  Rabi  oseillation  between  the  two 
states.  If  we  express  the  Rabi  frequency  as  ao  =  |  Ao|e"P,  the  solution  simplifies  to 

(6’o(t))  =  cos  a  con  11  sin  o  sin 

(C’J(t))  =  casino  cos  |Q0|(  +  it  *°  cosa sin  |Ji<>|<. 


the  dynamics  of  this  two-level  atomie  system  can  be  described  as  a  qubit  on  a  Bloeh 
sphere,  where  the  two  states  arc  represented  by  the  south  and  north  pole,  respectively 
(Fig.  24b).  The  initial  state  is  a  spin  coherent  state  on  the  Bloeh  sphere  characterized  by 
the  polar  angle  2( and  azimuthal  angle  ®.  The  evolution  of  the  atomic  state  interacting 
with  the  optical  field  over  time  T  corresponds  to  a  rotation  around  an  axis  n  in  the  x  -  y 
plane  by  the  amount  2 1 fio |  T'.  The  rotational  axis  h  passes  through  the  origin  and  makes 
an  angle  cp  with  the  x-axis. 

The  solution  found  above  is  formally  identical  to  the  ease  when  the  system  is  driven 


by  two  coherent  states,  but  with  two  critical  differences.  First,  the  Rabi  frequency  and  AC 
Stark  shifts  arc  given  by  the  ensemble-averaged  value  over  the  partially  coherent  field 
modes  whereas  such  averaging  is  not  necessary  for  pure  coherent  states.  Second,  the  va¬ 
riance  of  the  Rabi  frequency  is  non-zero  and  the  atomic  system  evolution  features  fluctu¬ 
ations  whereas  such  fluctuation  can  be  suppressed  for  the  coherent  state-driven  case.  This 
fluctuation  contributes  to  the  SNR  of  the  detection  scheme. 

Direct  measurement  of  the  atomic  population  after  the  interaction  provides  a  meas¬ 
ure  of  the  observable  |rj 2 ( r)|  rather  than  the  mutual  coherence  rj2(r).  In  order  to  measure 
rj2(r),  we  take  the  approach  of  adding  the  two  coherent-state  control  beams  to  the  system 
to  counteract  the  state  evolution  caused  by  the  two  optical  fields  under  investigation  (Fig. 
23b).  The  system  evolution  is  described  by  including  all  the  terms  in  the  Hamiltonian. 

The  resulting  solution  is  identical  in  form  as  above,  with  the  AC  Stark  shift  terms  and  the 
Rabi  frequency  replaced  by 

-  <$;  =  < 5.  +  <5*.  -  nr  =  fl0  +  n,. 

=  MViVal. 

fli  =  .^12«5o,r-  "V At, 

where  ao  and  a i  are  the  coherent  amplitudes  of  the  two  coherent  fields  and  ip  is  the  phase 
difference  between  them.  In  this  approach,  the  atomic  system  prepared  in  a  set  of  well- 
defined  states  (c.g.,  a  few  states  in  the  equator  of  the  Bloch  sphere)  interacts  with  the  two 
unknown  fields  and  the  two  control  fields  for  a  prescribed  time  T,  and  then  the  atomic 
population  is  measured  in  the  two  basis  states  |0)  and  |1)  with  close  to  unity  quantum 
efficiency.  The  amplitudes  and  relative  phase  difference  between  the  two  coherent  con¬ 
trol  fields  are  adjusted  as  the  experiment  is  repeated,  until  the  atomic  system  evolution  is 
completely  suppressed  (ftr=  0).  Under  this  condition,  the  coherence  of  the  control  fields 
infers  the  mutual  coherence  to  be  measured  (fto  =  -fii). 

The  sensitivity  of  this  approach  to  coherence  detection  is  limited  by  four  major  noise 
processes:  the  fluctuation  in  the  Rabi  frequency  due  to  the  fluctuation  of  the  optical  fields 
(“photon  noise”),  the  noise  associated  with  determining  the  coherence  of  the  control 
beams  (“optical  readout  noise”),  the  spontaneous  emission  of  the  atoms  during  the  Ra¬ 
man  process  (“spontaneous  noise”),  and  the  statistical  noise  from  detection  of  atomic 
states,  which  will  be  projected  to  one  of  the  two  basis  states  (“atom  shot  noise”).  The 
photon  noise  in  atomic  evolution  is  contributed  only  by  the  partially  coherent  fields  (and 
not  by  coherent  control  fields),  and  the  variance  averaged  over  the  measurement  time  T  is 
calculated  by  considering 

i  1  fT  (T  ,  2ir  n+1 

=_  /  <*i/  rf<a<n(<i)n(ia))2ifi2(i  +  _A  ). 

T  Jo  h\  T Au;  n 

The  variance  in  measurement  of  the  mutual  coherence  due  to  photon  noise  is  given 
by  cr  =|lj2|2(2^-/7’A<w)(»-i-l)//7.  The  shot  noise  corresponding  to  determining  the  mutual 
coherence  of  the  control  beams  is  given  by  21L  (ha>L  ITA),  where  7/  is  the  intensity  of  the 

laser  beam  interacting  with  the  atoms.  Under  the  balance  condition  It/Ai  ~  rnJ^/nA,  the 
optical  readout  noise  for  the  mutual  coherence  is  given  by 
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where  we  assume  a>\,2  -  coc  -  col.  For  the  spontaneous  noise,  we  require  that  the  sponta¬ 
neous  emission  rate  of  the  upper  state  of  the  atom  be  small  compared  to  the  Rabi  fre¬ 
quency  due  to  the  optical  fields,  so  that  the  “error”  from  spontaneous  emission  during  the 
coherence  measurement  is  small.  The  corresponding  noise  is  given  by 
<jsp  s(4;r2hco(  A/3/t/F).  where  A  -  2tt/ojc  is  the  central  wavelength  of  the  optical  fields. 

For  the  atom  shot  noise,  the  precision  with  which  one  can  determine  the  angle  of 
system  rotation  is  given  by  the  shot  noise  in  the  number  of  atoms  detected  in  each  state. 

The  corresponding  variance  in  mutual  coherence  is  given  by  of  athW'MiV?  /n.. 

where  <J,2  is  the  optical  transition  dipole  moment  between  the  state  |»)  and  |2)  .  The  over¬ 
all  variance  of  mutual  coherence  detection  in  the  quantum  coherence  detector  is  given  by 
<rQ=<rp+oi+o'p  +  cra. 


In  the  Gaussian  beam  geometry  with  a  beam  waist  of  Wo,  the  sensitivity  enhance¬ 
ment  factor  Qm.  m  of  /of  defined  as  the  improvement  over  classical  interferometric  ap¬ 


proaches 


Qnu 
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2  \J?)  \2(n+l)At  9AM»  +  1)Aw* 

+  (— _ _ )  . 
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where  rtr  =  nAcoT/2n  is  the  total  number  of  photons  detected  per  measurement  time  in¬ 
terval  T,  and  an  efficient  irradiancc  detector  was  assumed  (rj  =  1).  Figure  25  shows  Q„K  as 
a  function  of  nr  assuming  small  volume  cavity  with  high  finesse,  using  typical  values  of 

(J02  -  d\2  —  1.6x1 0 ~29Cm,  Aa>  ~  10'7coc,  A  =  780nm,  No  ~  lO'^m  3  is  the  density  of  the 
atoms,  Wo  =  1  pm,  and  J7  =  1.6*  10*. 


•6  4  2  0  2468 

log10nr 

Figure  25:  Enhancement  factor  Qmt  of  quantum  coherence  sensor  over  classical  interference  deteetors  plot¬ 
ted  as  a  function  of  nr-  Solid,  dashed  and  dotted  lines  represent  the  average  number  of  photons  per  mode  n 
=  10,  1  and  0.01,  respectively,  while  the  green,  blue,  red  and  black  curves  correspond  to  ^  values  of  1,  0.1, 

0.01  and  0.001,  respectively. 

At  low  nT,  the  sensitivity  of  the  quantum  coherence  deteetor  is  limited  by  the 
atom  shot  noise,  and  is  below  that  of  elassieal  interferometrie  deteetor.  As  the  average 
photon  number  per  mode  n  increases,  the  thermal  noise  in  the  optical  fields  increase  and 
the  quantum  eoherenee  deteetor  performance  improves  as  the  SNR  in  elassieal  measure¬ 
ment  degrades  (dotted,  dashed  and  solid  curv  es).  The  noise  performance  of  the  quantum 
deteetor  improves  as  nT  increases,  until  it  is  limited  by  the  optieal  readout  noise  or  the 
photonie  noise  limit.  In  this  limit,  the  dcpcndcnec  of  noise  on  nT  is  identical  for  the  quan¬ 
tum  and  elassieal  ease,  so  the  enhancement  faetor  Qme  remains  constant  as  a  function  of 
nT.  In  classical  interferometry,  the  noise  is  determined  by  the  photon  noise  corresponding 
to  the  irradianee  of  eaeh  optieal  field  (  m  i  and  ^2),  while  in  the  quantum  ease  the  noise 
contribution  only  comes  from  the  amount  of  mutual  coherence  present  between  the  two 
fields.  When  the  mutual  coherence  is  very  small  between  two  fields  with  low  average 
photon  number  per  mode,  the  quantum  eoherenee  deteetor  sensitivity  is  higher  by  a  faetor 

of  1  £ .  At  higher  values  of  nT,  longer  integration  time  increases  the  probability  of  spon¬ 
taneous  emission  and  the  sensitivity  is  determined  by  the  spontaneous  emission  noise. 

Through  our  theoretical  analysis,  we  have  effectively  proposed  a  “coherence  in¬ 
terferometer”  where  the  mutual  coherence  between  two  optical  fields  interferes  with  the 
mutual  eoherenee  of  two  eontrol  fields  (prepared  in  coherent  states).  The  mutual  cohe¬ 
rence  of  the  eontrol  fields  ean  be  manipulated  at  will,  and  we  balanee  the  eoherenee  inter¬ 
ferometer  so  that  the  mutual  coherences  of  the  pairs  of  fields  destructively  interfere.  This 
provides  a  means  to  quantitatively  estimate  the  mutual  eoherenee  in  the  signal  fields. 
There  has  been  some  debate  on  whether  the  noise  in  the  signal  fields  impacts  the  va¬ 
riances  of  coherent  rotations  of  the  atomic  system  that  accommodates  the  interference. 
Our  analysis  shows  that  the  impact  of  the  intensity  noise  from  individual  fields  do  not 
impaet  the  coherent  rotations,  but  it  was  pointed  out  that  the  variance  in  the  rotation  angle 
of  the  atomie  system  should  be  evaluated  by  considering  a  fourth-order  eoherenee  func- 


tion  of  the  Rabi  rotation  angle  Q  (which  will  be  an  8-th  order  coherence  function  of  the 
optical  fields).  Whether  such  higher  order  coherence  functions  impact  the  performance  of 
coherence  interferometer  has  not  been  analyzed  at  this  point. 

The  experimental  condition  that  has  to  be  satisfied  to  realize  such  coherence  inter¬ 
ferometer  seems  to  be  extremely  challenging,  mostly  because  the  coherent  interaction 
between  the  natural  fields  (the  signal  fields)  and  the  (far-detuned)  atomic  resonance  is 
necessarily  very  weak.  Spontaneous  emission  will  lead  to  deeoherenee  of  the  coherence 
interferometer  unless  the  interaction  can  be  dramatically  enhanced.  The  use  of  high  fi¬ 
nesse  optical  cavity  could  address  this  issue,  but  w  ill  lead  to  constraint  on  the  bandwidth 
of  the  signal  fields.  One  has  to  distinguish  the  effect  of  severe  filtering  from  the  cavity 
(that  enhances  coherence)  and  the  improved  sensitivity  of  our  coherence  detector  due  to 
enhanced  interaction  by  the  cavity.  Our  assumptions  operate  in  the  regime  where  the  cav¬ 
ity  bandwidth  is  still  large  compared  to  the  bandwidth  of  the  signal  fields  and  demon¬ 
strate  the  benefits  of  the  mutual  coherence  interferometer.  The  remaining  conditions  of 
the  experiments,  including  the  atomic  density  and  coherence  times  required,  makes  this 
an  extremely  challenging  experiment  to  perform  in  reality. 

Compressive  Holography 

Compressive  measurement  was  the  second  emphasis  of  the  coherence  measure¬ 
ment  portion  of  this  program.  We  realized  in  early  2009  that  compressive  inference  strat¬ 
egics  developed  in  support  of  CASSI  imaging  could  be  applied  to  form  tomographic  im¬ 
ages  from  single  frame  Gabor  holograms.  We  demonstrated  this  capacity  in  stand-off 
lcnsless  imaging  in  [3].  As  a  second  demonstration,  we  demonstrated  3D  reconstruction 
of  a  microscopic  object  from  a  single  2D  Gabor  hologram  produced  by  the  microscope 
with  a  point  source.  For  compactness,  these  microscopes  typically  use  the  Gabor  geome¬ 
try  rather  than  the  Leith-Upatneik’s  geometry.  This  causes  the  undesirable  autocorrela¬ 
tion  and  twin-image  problems.  Hence,  our  goals  are  two-fold:  1)  to  solve  the  underdeter- 
mined  problem  of  obtaining  a  3D  tomographic  reconstruction  from  a  2D  hologram,  and 
2)  to  remove  the  autocorrelation  and  twin-image  problems  that  often  severely  deteriorate 
the  visual  quality  of  the  reconstruction. 

Figure  26  shows  the  system  schematic  of  compressive  holographic  microscope. 

A  coherent  light  is  spatially  filtered  and  collimated  to  make  a  plane  wave  and  this  plane 
wave  is  focused  to  make  a  point  source  by  a  microscope.  From  this  point  source,  a  sam¬ 
ple  is  illuminated  and  magnified  at  the  FPA.  The  recorded  interferogram  on  FPA  is  the 
Gabor  hologram  with  the  spherical  wave. 


Fig.  26.  Schematic  of  compressive  holographic  lcnsless  microscope 


In  the  Gabor  geometry,  the  intensity  measurement  on  FPA  is  expressed  as 

=  |/l(.t,.v,z)|'  +  |/'.(v. y; y; z)l-(\,y\z)  +  a(x,  v;z)e’(x, y;z). 


(1) 


where  A  is  the  reference  Field  with  spherieal  phase  profile  from  the  point  souree  at  r=z0 
and  E  is  the  scattered  objeet  field.  This  scattered  field  from  a  3D  object  density//  ean  be 
defined  by  the  first  Bom  approximation  as 


(2) 


/'.’(v,r;  r )  =  JjJ  dx'dy'dz' p{x' ,  y\z']A(x',y';z') i(x  -  x',  y 


Here,  h  is  a  propagation  kernel  and  Fourier  transform  of  this  propagation  kernel  is  given 
by 


(3) 


The  deteetor  plane  is  positioned  at  z= 0.  Sinee  this  propagation  is  linear  operation,  the  Ga¬ 
bor  hologram  ean  be  represented  as 


(4) 


g  =  2  Re{///" }  *  e  +  n , 


where /is  the  discrete  form  of  /?.  e  denotes  the  autocorrelation  of  A  and  E,  and  n  is  addi¬ 
tive  noise.  In  Eq.  (4),  by  enforcing  a  sparsity  constraint  on  the  total  variation,  the  discrete 
objeet  density/is  estimated  as 


/  =  argmin(j|g-///||‘+||/|rJ. 


(5) 


Here,  ||/fln  is  total  variation  norm  of  veetor / and  defined  as 

mu  =ZZZI^>J’ 


(6) 


z  x  y 


where  f.  denotes  a  2D  plane  of  the  3D  objeet  dataeube. 

We  demonstrated  compressive  holographic  microscopy  using  a  He-Nc  laser  with 
632. 8nm  wavelength  as  a  light  souree.  A  mieroseope  with  0.65NA  is  applied.  The  focal 
plane  array  is  a  Lumenera  CMOS  sensor  with  I280xi040 resolution,  5.2^1  pixel  pitch  and 
10  bit  digitization.  Figure  27  shows  the  experimental  setup  and  a  structure  of  the  lensless 
mieroseopy.  The  positions  of  the  mieroseope  and  the  sample  are  adjustable  with  3-axis 
stages. 

In  this  system,  the  magnification  is  determined  by  the  ratio  of  effective  distances 
from  the  point  source  to  the  sample  and  from  the  point  source  to  the  FPA.  The  numerical 
aperture  (NA)  of  this  system  is  a  function  of  pixel  pitch  of  FPA,  the  wavelength,  and  the 
magnification.  As  the  most  simple  magnification  system,  the  NA  of  4 / optics  is  given  as 


(7) 


where  M  is  magnification,  and  p  is  the  pixel  pitch  of  FPA.  The  NA  of  this  system  is  larg¬ 
er  than  that  in  a  magnification  system  with  4f  opties  even  though  the  wave  recorded  on 
FPA  has  a  spherieal  phase.  The  optical  magnification  is  set  as  4.5  and  then  the  NA  is 
0.27,  whieh  yields  axial  resolution  a_.  =2 a/na2  =16.9/™ and  transverse  resolution 

A,  =  X/2NA  =!.15/»n. 


As  a  sample,  a  water  flea  with  a  1 ,5mm  size  is  measured  using  a  snap  shot.  Figure  28 
shows  a  raw  Gabor  hologram  and  reconstruction  images.  The  transverse  slices  at  various 
ranges  with  baek  propagated  fields  and  eompressive  holography  are  shown  respectively 


in  Figs.  28(b)  and  (c),  where  each  interval  between  adjacent  slices  is  0.3 mm .  In  the  slices 
with  compressive  holography,  the  antenna  structures  are  clearly  shown  with  the  efficient 
removal  of  the  autocorrelation  and  twin-image  terms  in  comparison  of  the  images  with 
back  propagated  fields. 


(a) 


(b) 


Figure  27.  (a)  Photograph  of  experimental  setup  and  (b)  structure  of  holographic  microscopy 

We  demonstrated  the  feasibility  in  obtaining  a  3D  reconstruction  of  a  microscopic 
object  from  a  single  2D  Gabor  hologram  produced  by  a  holographic  microscope.  In  Figs. 
28(b)  and  (e),  the  reconstruction  images  with  back  propagated  fields  shows  severely  dete¬ 
riorated  images,  but  the  images  w  ith  compressive  holography  shows  the  clear  antenna 
structure,  since  the  linear  twin  image  term  and  splits  the  non-linear  autocorrelation  term 
in  the  measurement  plane  are  significantly  removed.  However,  there  arc  only  small  dif¬ 
ferences  among  slices  of  dataeube  even  though  the  distance  between  slices  arc  enough 
longer  than  axial  resolution.  We  conjecture  that  the  reason  is  that  the  object  is  rather  thick 
in  comparing  the  transverse  resolution  and  the  radiation  scattered  from  object  has  narrow 
distribution.  The  edge  in  the  body  of  the  specimen  is  only  observable  in  reconstruction 
images,  since  the  Gabor  hologram  is  hard  to  measure  the  part  of  object  with  uniformly 
low  transmission  in  large  area.  The  reconstruction  errors  and  the  effects  of  noise  usually 
results  from  the  insufficient  number  of  measurements.  It  w  ould  be  desirable  to  combine 
this  system  with  synthetic  aperture  technique  or  structured  illumination  for  increasing  the 
number  of  measurements. 


Fig.  28.  (a)  Raw  Gabor  hologram  measuring  a  water  flea,  (b)  transverse  slices  at  various  ranges  with  back 
propagated  fields,  and  (c)  transverse  slices  at  various  ranges  with  compressive  holography. 
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Sc  Hoon  Lim 
Kcrkil  Choi 
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15,  14013-14027  (2007) 

5.  Fernandez,  C.A.,  et  al..  Fluorescence  microscopy  with  a  coded  aperture  snapshot 

spectral  imager.  Proc  SPIE,  2009.  7184-33:  p.  1-11 

6.  Christy  Fernandez  Cull,  Kerkil  Choi,  David  J.  Brady,  and  Tim  Oliver,  "Identifi 

cation  of  fluorescent  beads  using  a  coded  aperture  snapshot  spectral  imager," 
Appl.  Opt.  49,  B59-B70  (2010) 

7.  Cull,  C.F.,  et  al.,  Direct  fluorescence  spectral  feature  identification  using  a 

coded  aperture  snapshot  spectral  imager.  Optics  Express,  2009. 

(S’.  Brown,  B.,  et  al.,  Growth  of  aligned  bamboo-structured  carbon  nanotubes  via 
platinum  catalyst,  in  preparation  for  Nanoletters. 

7.  Interactions/Transitions: 

a.  Participation/presentations  at  meetings,  conferences,  seminars,  etc. 

K.  Choi,  R.  Horisaki,  D.  Marks,  and  D.  Brady,  "Coding  and  Signal  Inference  in  Com¬ 
pressive  Holography,"  in  Computational  Optical  Sensing  and  Imaging,  OSA  Technical 
Digest  (CD)  (Optical  Society  of  America,  2009),  paper  CThA5. 

S.  Lim,  R.  Horisaki,  K.  Choi,  D.  Marks,  and  D.  Brady,  "Experimental  Demonstrations  of 
Compressive  Holography,"  in  Computational  Optical  Sensing  and  Imaging,  OSA  Tech¬ 
nical  Digest  (CD)  (Optical  Society  of  America,  2009),  paper  CThA6. 

A.  Wagadarikar,  D.  Marks,  K.  Choi,  and  D.  Brady,  "Compressive  Coherence  Sensing,"  in 
Computational  Optical  Sensing  and  Imaging,  OSA  Technical  Digest  (CD)  (Optical  So¬ 
ciety  of  America,  2009),  paper  CTuA6. 

C.  A.  Fernandez,  S.  Lim,  B.  D.  Guenther,  D.  J.  Brady,  and  S.  T.  McCain,  "High- 
Throughput,  Multiplex  Aperture-Coded  Raman  Spectrometer  for  Biomedical  Diagnos¬ 
tics,"  in  Biomedical  Optics,  OSA  Technical  Digest  (CD)  (Optical  Society  of  America, 
2008),  paper  BTuF15. 

Materials  Research  Society,  Materials  Research  2008,  North  Carolina,  Organizing  Com¬ 
mittee,  November  18.  2008 

Bouehet  Society  (for  underrepresented  minority  graduate  students  in  the  natural  sciences 
engineering,  and  mathematics),  Duke  University,  Poster  Presentation,  Billyde  Brown  et 
al. 


Materials  Research  Society,  Materials  Research  2008,  North  Carolina,  Organizing  Com¬ 
mittee,  November  18.  2008,  Poster  Presentation,  Billyde  Brown  et  al. 

b.  Consultative  and  advisory  functions 

Panelist,  STaR  (Seienee,  Technology  and  Research)  Symposium,  April  14-15, 
Charleston,  West  Virginia 

c.  Technology  Assists,  Transitions,  and  Transfers. 

Coded  aperture  snapshot  spectral  imaging  systems  derived  from  research  under 
this  program  have  been  developed  in  prototype  form  for  several  users  by  Applied  Quan¬ 
tum  Technologies,  Ine.  A  competitive  test  of  this  technology  with  conventional  spectral 
imagers  is  scheduled  with  the  National  Geospatial  Agency  for  June  2010.  A  CASSI  sys¬ 
tem  is  also  scheduled  for  delivery  at  the  Peabody  Museum  at  Yale  University  in  2010. 
CASSI  has  proven  to  be  a  unique  and  useful  approach  for  rapid  spectral  image  capture, 
we  anticipate  continuing  transitions.  Compressive  holography  studies  under  this  program 
led  to  a  new  DARPA  MTO  seedling  program  in  holographic  remote  sensing  (collabora¬ 
tion  with  the  University  of  Rochester.)  Compressive  holography  results  were  also  pre¬ 
sented  at  the  industry  day  for  DARPA’s  Fine  Domain  Optical  Surveillance  (FDOS)  pro¬ 
gram.  Compressive  holography  is  a  component  of  system  design  for  several  FDOS 
projects. 

8.  New  discoveries,  inventions,  or  patent  disclosures. 

CNT  growth  using  Pt  catalyst  has  crossover  application  in  the  area  of  nerve  sti¬ 
mulation.  Thus,  an  Invention  Disclosure  has  been  submitted  entitled:  “ Biocompatible 
Nanostruetured  Coating  To  Reduce  Risks  Of  Physiological  Harm  And  Material  Damage 
Of  Implantable  Electrodes  During  Chronic  Nerve  Stimulation  Therapy” 

9.  Honors/ Awards: 

Nan  Zheng  passed  PhD  preliminary  exam  4/21/08 

C.  Fernandez  (National  Defense  and  Engineering  Graduate  Fellowship) 

A.  Wagadarikar  (Natural  Sciences  and  Engineering  Research  Council  [NSERC]  of  Cana¬ 
da  Fellowship) 
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